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Abstract 

A generalized Fokker-Planck equation is derived to describe particle kinetics in specific situations 
when the probability transition function (PTF) has a long tail in momentum space. The equation is 
valid for an arbitrary value of the transferred in a collision act momentum and for the arbitrary mass 
ratio of the interacting particles. On the basis of the generalized Fokker-Planck equation anomalous 
diffusion in velocity space is considered for hard sphere model of particle interactions, Coulomb 
collisions and interactions typical for dusty plasmas. The example of dusty plasma interaction is 
peculiar in way that it leads to a new term in the obtained Fokker-Planck-iike equation due to the 
dependence of the differential cross-section on the relative velocity. The theory is also applied to 
diffusion of heavy particles in the ambience of light particles with a prescribed power-type velocity 
distribution function. In general, the theory is applicable to consideration of anomalous diffusion in 
velocity space if the typical velocity of one sort of particles undergoing diffusion is small compared 
to the typical velocity of the background particles. 

PACS number(s): 52.27.Lw, 52.20.Hv, 05.40.-a, 05.40.Fb 

Introduction 

Interest to anomalous diffusion in coordinate space is explained by its great variety of applications, 
i.e., semiconductors, polymers, some granular systems, plasmas under specific conditions, various 
objects in biological systems, physical-chemical systems, and others. 

Non-linear time dependence < r 2 (t) > of the mean-squared displacement has been experimentally 
observed, in particular, under essentially non-equilibrium conditions and in some disordered systems 
[TJ [2]. As is well known there are two types of anomalous diffusion in coordinate space which 
are referred to as superdiffusion < r 2 (t) >^ t a (a > 1) and subdiffusion (a < 1) |3J. For 
description of these two anomalous diffusion regimes a number of efficient models have been 
proposed. The continuous time random walk (CTRW) model by Scher and Montroll |1] for the 
process of subdiffusion is a basis for understanding photoconductivity in strongly disordered and 
glassy semiconductors. The Levy-flight superdiffusion model [5] describes such phenomena as self- 
diffusion in micelle systems [3] as well as reaction and transport in polymer systems [8]. Another 
application of the model is studying stochastic behavior of financial market indices [5] . The method 
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of so-called fractional differential equations in coordinate and time spaces has been successfully 
developed to treat both cases of anomalous diffusion [9] . 

Recently a more general approach has been proposed in [10J which not only reproduces the 
results of the standard fractional differentiation method (whenever it is applicable) but also allows 
to tackle more complicated cases of anomalous diffusion. In [11], for example, this approach is used 
for studying diffusion in a time-dependent external field. 

Some aspects of anomalous diffusion in velocity space have been considered in papers |13|-|18]. 
However, compared to anomalous diffusion in coordinate space anomalous diffusion in velocity space 
is still poorly studied but it is actively attracting more and more attention [19|-[21j. 

In this paper the problem of anomalous diffusion in momentum (velocity) space is considered 
in detail in the spirit of the approach proposed in [10] for the diffusion in coordinate space. Extra 
terms (cp. [Hj]-[2I]) of the generalized Fokker-Planck equation are found. The results are applied to 
hard sphere model of collision, Coulomb systems and dusty plasma. In dusty plasma the additional 
terms are numerically small but not negligible due to to velocity dependence of the cross-section. 

The paper is organized as follows. Diffusion in velocity space for the cases of normal and 
anomalous behavior of the PTF is presented in Section 1. Starting from the Boltzmann-type PTF 
we arrive at explicit expansion PTF into a series which is applicable not only to the Boltzmann-type 
processes but also to a wide class of processes with other types of PT-functions. The generalized 
Fokker-Planck equation is presented in Section 2. Special cases of anomalous diffusion with the 
specific power- type distribution function of ambient light particles are analysed in Section 3. In this 
section the cases of hard sphere and Coulomb models of diffusion are considered. Subsection 3.3 is 
devoted to the example of anomalous diffusion in dusty plasma where the differential cross-section 
depends on velocity of the colliding particles. In this case the additional term in the generalized 
kinetic equation is not negligible. 



1 Diffusion in velocity space on the basis of the master-type 
equation 

Let us consider diffusion in velocity space using an approach based on the corresponding master 
equation for the distribution function /(p, i) which describes the balance of particles coming at and 
from the point p at the instant t. The master equation written in velocity space reads (see, e.g., 



T =Jdq {W(q, p + q)/(p + q, t) 

-W(q,p)f(p,t)}. (1) 

The probability transition function VF(q, p') describes the probability that a particle with momentum 
p' passes from the point p' in the velocity phase space to the point p per unit time by transferring 
momentum q = p' — p to the surrounding medium. Under the assumption that the characteristic 
transfered momentum q is much smaller than p equation (fTJ) can be expanded into a series with 
respect to q to the second order. The expansion leads to the usual form of the Fokker-Planck 
equation for the density distribution function /(p,i) 



2 



4tlzH-J>-[A a (p)f(p,t) 



dt 



■^(B af) (p)f(p,t)) 



(2) 



where 



A a{p) = J q a W(q,p)dq, 
B a p(p) = \ J QaQpW(q,p)dq. 



(3) 
(4) 



The coefficients A a and B a p are responsible for the processes of friction and diffusion, respectively. 
The subscripts a and /3 correspond to the coordinate axes. 

For the Boltzmann elastic collisions the PTF-function W(q, p) has been found in |23[ and 
can be represented in the form [TU] 



W(p,q) = — 
A* 



duS uq 



da 
do 



arccos (1 — 



2/i 2 u 2 



),u 



2(i. 

A(u + v). 



(5) 



Here p = Mv is the momentum of a scattered particle, \i = Mm/(m + M) is reduced mass 
and ft is the distribution function for the background scattering particles. The relative collision 
velocity is denoted with u and da(x,u)/do is the scattering differential cross-section calculated in 
the coordinate system in which a scattering particles with mass m is fixed throughout the whole 
act of collision [25] [26] . The scattering angle \ = arccos[l — q 2 / (2[i 2 u 2 )] stands for the rotation 
of the relative velocity u. Assuming that the vector q is directed along the axis z and performing 
integration with respect to the projection of u on the axis z one arrives at the more explicit form 
of the equation for the PTF-function 



w(q.p) = i/ du ^ M« 2 =«2+«i 
, [, ^2 2 q ,v (q ,v ) 2 

L A* 1 



(6) 



where u z = —q/2fi , and mean projections of the vectors u and v, respectively, on the plane 
xy perpendicular to the axis z. 

For the scattering cross-sections independent on u (e.g., for hard sphere or Coulomb interactions) 
equation ([6]) is reduced to the following equation by changing the variable of integration u = u^+v^ 



7r da 

W{q.,p) = — — 
[i £ q do 



fb 



u 2 



q 



q-y 



(q-v) 2 



(7) 



du 2 



Now let us consider the case when the characteristic velocity of particles diffusing in velocity 
space is small compared to the characteristic velocities in the surrounding medium (e.g., [22 ). 
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For the Boltzmann-type collisions it means that there are two types of particles - heavy particles 
undergoing diffusion with the characteristic velocity v ~ p/M and light particles characterized 
by velocity u and u > «. In the zero order approximation with respect to p when we omit any 
dependence of the PTF on p (i.e., W(q, p) = W(q, p = 0) = W{q)) the coefficients A a and B a p 
yield the equations 



A a = 0, 



S f 

B a p = S aP B = ~f- q 2 W(q)dq 7 



(8) 



where 8 a p is the Kronecker delta. The neglect of p-dependence is not correct. Usually the coefficient 
A a for the Fokker-Planck equation is found assuming that the stationary distribution function is 
Maxwellian. This assumption leads to the standard relation for the coefficients MTA a (p) = p a B 
where M is mass of the diffusing particle and T is temperature of the particles in equilibrium. The 
relation is analogous to the Einstein relation in coordinate space. 

For the systems far from equilibrium, for example, for slowly decreasing PTF, the approach 
above is not applicable. In this case rigorous calculation of the coefficients A a and B a p requires 
the next order of approximation with respect to the smallness of p [T5]. As the function W(q, p) 
is a scalar it can be regarded as a function depending only on variables q, q-p, p 2 . We imply that 
the PTF is such a function of q • p and p 2 that it can be analytically expanded into a series in the 
vicinity of the point q • p = p 2 = 0. Note that in the expansion of the PTF we leave only those 
terms which have p to the power no higher than 2. Expansion of W(q, p) into a series to the second 
order is given by 

W(d,p)~W(q)+W'(q)(ci-p) 

1 - (9) 
+ -W"(q)(c l -p) 2 + W'(q)p 2 



where 



W'(q) -- 

W"(q) 

W'(q) 



dW(q,q-p,p 2 ) 



9(q-p) 
d 2 W(q, q-p,p 2 



<9(q-p) 2 
dW(q,q-p,p 2 ) 



dp 2 



q-p— 0, p 2 =0 

(10) 

q-p=0,p 2 =0 
qp=0,p 2 =0 



Basically, the developed approximation is true if the typical velocity of one sort of particles undergoing 
diffusion is small compared to the typical velocity of the background particles. Let us note that for 
the differential cross-sections independent on the relative velocity u the coefficient W'(q) = as it 
can be seen from equation ©• In general, this term is not negligible and it should be taken into 
account (it was omitted in [Hfl I21| ). 
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Then with the necessary accuracy we find the coefficient A a 

4q(p) = J qaqpPpW' (q)dq 

= Pa J q a q a W'{q)dq (11) 
= y / q 2 W\q)dc l . 

Remarkably, if the equality W'(q) — mW (q) / (2M fiT) is fulfilled for the function W(q, p) the 
following relation for the coefficients A a and B takes place 

—MTA a (p) = p a B (12) 

TO 

which turns into the usual Einstein relation in the limit m <C M 



MTA a (p)= Pa B. (13) 

Straightforward differentiation of W(p, q) in Eq. (p^J) with respect to the variable pq which is 
assumed independent in the case of the equilibrium Maxwellian distribution leads to equality (|12| 
and we arrive at the ordinary Fokker-Planck equation in velocity space with the constant diffusion 
coefficient D = B/M 2 and the constant friction coefficient (3 = mB/ ' \iT = D(m + M)/T which 
satisfy the Einstein relation in the limit to <C M. 

For systems close to equilibrium the PT-function is calculated and the appropriate Fokker- 
Planck equation is discussed in detail in review |23| . In case of hard-sphere interactions and the 
equilibrium Maxwellian distribution /& equation §5§ is reduced to the equation which is consistent 
with the result given in papers |23[ 124] . However, even for quasi-equilibrium regimes when long 
tails of the PTF are absent the consideration in [53] is restricted to hard sphere interactions. Non- 
equilibrium forms of the PTF and changes in the structure of the Fokker-Planck equation for such 
systems are considered in [19, 20J. In the following section this consideration is extended. It brings up 
the importance of of the additional term in the expansion of the PTF with respect to the variable p 2 
(see Eq. ©). The chapter also gives correct coefficients in the generalised Fokker-Planck equation. 



2 Generalized Fokker-Planck equation 

In some non-equilibrium (stationary or non-stationary) systems PTF can have a long tail as a 
function of q. For such systems the ordinary Fokker-Planck equation @ is not valid as the kinetic 
coefficients diverge at large values of q. 

Now let us refer to the expansion of the PTF given by equation ([9]) making no assumptions 
as for the dependence of the differential cross-section on the relative velocity u. It results in an 
additional (cp. [T{?1 120|) term proportional to the first derivative of the PTF with respect to p 2 . 
Substituting this expansion into Eq. ([T]) and making use of the Fourier transformation we arrive at 
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the following equation 



dm 

dt 



= J W(q) [/(p + q) - /(p)] ^dpdq 



+ / W"(g)[(q-(p + q))/(p + q) 



-(q-p)/(p) 



e *p.s 



dprfq 



(2tt) 3 

+ / ^"(<?) [(q-(p + q)) 2 /(p + q) 



(14) 
e ip-s 

-(q-p) 2 /(p)] (^j3 d P rf q 

+ /" W"(g) [(p + q) 2 /(p + q) 
e ip-s 

where /(s) is the Fourier image of the function /(p). We omitted the arguments t and q of the 
distribution function to make the notation brief. Introducing the coefficients 



A(s) = 
B a (s) 



e -i(q-s) _ 1 
i 



(q-s) 



C a0 {s) = -\ J 



W(q)dq, 



-i(q-s) 



j(q-s) _ 2 



W'(q)dd 
W"(q)dci, 



25(s) 



e -i(q-s) _ X 



W'(g)dq 



we come to the Fokker-Planck-like equation 

df(s 



dt 



A(s)/(s) + B a (s) 



g/(g) 
<9s Q 



+ C a/3 (s) 



After a number of calculations it can be shown that 



£(s)A/(s). 



A(s) = 



= 4tt / 
Jo 



sin(gs) 



- 1 



W(«)dg, 



(15) 



(16) 



(17) 



B Q (s) = s a B(s), 
, x 4tt / >00 2 



cos(qs) 



sin(gs) 



qs 



W'(q)dq, 



(18) 
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C a p(s) = s a spCi{s) + s 2 5 a pC 2 {s), 



CM - ^ 



sin(gs) cos(qs) 

3gs 
sin(gs) 



q 2 s 2 



C 2 {s) 



-in 

~72 



q 3 s 3 
1 

3 " 



w"(g)dg, 

cos(gs) 



q 2 s 2 



sin(gs) 



g 3 s 3 



W"(«)dg, 



£(s) - s 2 C 3 (s), 



/ ^ 


' sin(gs)" 


/o 


gs 



W'(g)dt/. 



(19) 
(20) 

(21) 

(22) 
(23) 



For an isotropic function /(s) = f(s) the Fokker-Planck-like equation (fT6| is reduced to equation 

df(s) 



dt 



A(s)f(s) 



+ {B(s) + 2C 2 (s) + 2C 3 (s)]s 



df(s) 



(24) 



[Ci(s)+C 2 (a) + C 3 00] 



9s 2 ' 



This equation is a generalisation of the usual Fokker-Planck equation which implies certain smallness 
of the transferred momentum q. The advantage of our equation is that it is valid for an arbitrary 
value of q. This virtue lets us go far beyond the scope of the phenomena described by the ordinary 
Fokker-Planck equation. 

However, it is instructive to match Eq. and the usual Fokker-Planck equation. The latter 
tackles the problem of diffusion of heavy particle in a gas of light particles. Basically, Eq. (I24p 
should be reduced to the Fokker-Planck equation in the limit m/M — ► and for small values of q. 
First of all, let us note that in the case when the PTF and the functions W(q), W'(q), W'(q) and 
W"(q) strongly decrease for large values of q the exponents in the integrals in the functions A(s), 
B(s), Ci(s), C 2 (s) and 63(5) can be expanded into a series as follows 





s 2 
6 j 


f q 2 W(q)dq, 


B(s)~ 


~\I 


q 2 W'(q)dq, 


C X ( S ) ~ 


1 

~i~5 


J q 4 W"(q)dq, 


C 2 (s) ~ 


1 

~30 


J q 4 W"(q)dq, 


C 3 (s) ~ 


II 


q 2 W'{q)dq. 



(25) 
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Then we can set the coefficients C± , C2 and C3 equal to zero as they are of next order of smallness 
compared to the coefficients A and B with respect to the small parameter m/M. The fact can be 
deduced from the expression for the PTF ([6]) . And the simplified kinetic equation in velocity space 
based on the PTF (which is non-equilibrium in the general case) yields 

where A a = -1/6 / q 2 W(q)dq. 

The stationary solution of this equation is given by 

(27) 



2B 



/(s) = Const • exp 
The respective normalized stationary momentum distribution reads 

(28) 



NB 3 / 2 r Bp 2 



/(p)= (2^4^ eXP 



3 Special cases of anomalous diffusion in velocity space 

Let us calculate the kinetic coefficients in the generalized Fokker-Planck equation for the special 
cases of anomalous diffusion. In our model diffusing heavy particles with mass M interact with light 
particles of the surrounding medium with mass m <C M and with power-type velocity distribution 

A(u + v)= 3 " b (29) 
u '(u + v)^ 

where uo means the characteristic velocity, rib is the constant which has the dimension of the density 
in the coordinate space but since the function /& can not be normalized we can not think of rib 
as of a density. Impossibility to normalize the distribution function /& is explained with the fact 
that the normalizing integral diverges either at small or at large values of the argument of In 
reality normalization can always be fulfilled. At smal values of the argument the power-law (1291) of 
the distribution function is violated. But it is highly unlikely to obtain an analytical expression for 
W(% p) for realistic distribution functions. 

3.1 Hard-sphere interactions 

Interactions in systems of hard-sphere particles of radius a are described by the differential cross- 
section da/do = a 2 /A. Substitution of the specified distribution function and cross-section into 
equation §5§ gives 



riba 2 



«(«. + 



For 7 > 1 the integral (|3"D|) converges and the result of integration reads 

, . Trn h a 2 (£- 1) 2 ~ 2 ^ , . 

(2^) 4 - 2 7( 7 - l)g 2 T-lu5 7 
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where £ = 2/i(q • v)/g 2 . Expanding equation (f3"Tj) into a series according to © we leave only the 
first three terms 

W (q) = ^ = *° (32) 



27+1 ' 



M(2/i) 3 ~27g2 7 + l u 3-27 q 

2M2 7 -2) 
i?i - ^ i?o, 

^"( g )= 27m 6 a 2 (2 7 -l) _ i? 2 



M2(2/i) 2 -27 g 2 7+ 3 u 3-27 g 2 7 +3 ' 

(2/x) 2 (2 7 -2)(2 7 -l) D 

^ " M 2 «0> 



(33) 



(34) 



VK'(g) = 0. (35) 

Now let us substitute VF(g) into equation (jTTJ) which introduces the coefficient A(s) and discuss 
the convergence criteria. 



A(s) = 47r 

= 4tt.Ro 



pOQ 

/ 92 


'sin(gs) " 


/o 


qs 



dq 
g 2 7-3 



W(<z)dq 

(36) 



sin (ga) _ 
9« 



The convergence of the integral in equation (|36l) is guaranteed when the inequality 2 < 7 < 3 
takes place. Whereas 7 < 3 gives the convergence for small values of q (q — > 0) the convergence for 
q — > 00 is provided at 7 > 2. 

The values of the parameter 7 at which the coefficient B(s) remains finite can be found directly 
from the definition of B(s) (Eq. (fT5)0 . It is clear to see that condition for convergence of the integral 
in equation for the coefficient B(s) for small values of q is 7 < 2 and for large values of q is 7 > 1/2. 
Finally, finiteness of the coefficients Ci(s) and 62(3) follows from the inequality 7 < 2 for small 
values of q. For large values of q the coefficient C\ is finite at 7 > 1 and the coefficient C2 is finite 
at 7 > (Eq. pip!])). 

To sum up, in the case under consideration, the convergence of the coefficients A, B, and 
for large values of q is defined only by the convergence of the coefficient A, which means 7 > 2. 
For small values of q it is sufficient to guarantee convergence for the coefficient B(s) with 7 < 2. In 
other words, for purely power- type behavior of the function /&(£) the simultaneous convergence of 
the coefficients A, B and does not exist. However, in real physical models, the convergence of 
the coefficients at q —¥ can be achieved, for example, due to the non-power-type behavior of the 
PTF W for small values of q (compare with examples of anomalous diffusion in coordinate space 
|12j). Therefore, in the system of hard-sphere heavy particles undergoing anomalous diffusion in the 
medium of light particles with power-type distribution at large values of the variable q the kinetic 
coefficients exist at 7 > 2. This result approves of the qualitative conclusion in paper |19| but the 
numbers do not coincide as we started from another more explicit expression for the PTF ([5]). 



9 



3.2 Coulomb collisions 

By analogy, we can consider a system of particles characterized by Coulomb interactions. Formally, 
the Coulomb interactions are described by a differential cross-section 

da Coul _ ( Ze 2 \ 2 1 ^ 4Z 2 eV 
do \2trn 2 ) sin 4 X /2 g 4 ' 1 ' 

where Z is the charge number and the scattering angle \ = arccos(l — q 2 /2^,u 2 ) (see Eq. ([5])). 
From the mathematical point of view analysis of the kinetic coefficients in the case of Coulomb 
collisions is no more difficult than the case of hard-sphere collisions as the Coulomb cross-section 
is also independent on the relative velocity it. The only difference is the power in the dependence 
of the PTF on the transfer momentum q 

Performing integration we arrive at the equation 

, , 47rZ 2 e 4 n b (e - I) 2 " 27 , , 

Wcoui P, q = — 5-^7, 39 

y ' (2/i) 2 - 2 ''(7-l) (Z 2 7+3 u 3-27' V I 

Expanding the PTF into a series we obtain the functions W^Coui(<?), T^coui^) and WQ oul (q) 

„ T , , 47rZ 2 e 4 n b u 2 1 7 ~ 3 K , . 

Wc ° M = (2 M )'-^( 7 -l)^ +3 = { ) 

~ , _ 87rZ 2 e 4 n b Uo 7 " 3 K x 



= 2,(27-2) 
M 



- „ _ 8^ 2 e 4 n b (2 7 - l)^ 7 " 3 _ # 

^CoulW — 



M 2 (2^ 2 V 7+7 9 2 7+7 ' 
(2 M ) 2 (2 7 - 2)(2 7 -l) ^ 
^ 2 = M 2 A °' 

^Coui(9) - 0. (43) 

Reasoning by analogy with the previous section concerning hard-sphere interactions we can draw 
a conclusion that the coefficients A(s), B(s) and 61,2(3) are finite at large values of q provided the 
inequality 7 > is fulfilled. We assume that the convergence of the integrals in these coefficients 
at small values of q is due to various neglected in our model factors, such as screening effect, 
etc. The asymptotic behavior of the PTF and the respective derivatives of the PTF for large q is 
similar to one obtained in [21]. However, the coefficients are calculated on the basis of the explicit 
representation for the PTF ((6]). 
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3.3 Interactions typical for dusty plasmas 

A more complicated and interesting case is interaction of dusty particles with electrons and ions. 
The differential cross-section for this type of interaction depends on the relative velocity of the 
particles unlike the cases of Coulomb and hard-sphere models of collisions. This fact leads to the 
dependence of the PTF on p 2 . To demonstrate this let us refer to the the differential cross-section 
for grains-electrons interactions [57] . The respective differential cross-section yields 

2 / 2e 2 \ 2e 2 

TP I 1 j ' 2 < 1 ' 

a c (u) = { V / meUj PmcU (44) 
0, 2>1, 

a i (u) = Trp 2 (l + -^), (45) 

where m e (n is electron (ion) mass, p is radius of a grain, e is the elementary charge. The ions are 
assumed singly charged. 

Let us consider electron-grain collisions and rewrite Eq. (|44p in a more convenient way 

rr,.Un={ ^0 " (46) 



0, u < A, 

where A 2 = 2e 2 /pm e . Along with the previous special cases we assumed scattering particles 
(electrons) have a power- type velocity distribution (|29p . Substituting the announced cross-section 
and the distribution function into the PTF ((5]) we arrive at 



4^m 'q Ju>A \ u 
x 5 u z + — 1 



(47) 



2 AV ( u 



2- ' 



Analyzing the integral we can conclude that for large values of q/2p > A the inequality u > A 
is guaranteed. In our paper we restrict our consideration only with large values of q. The point is 
that the power-type distribution function is realistic only for large values of q. For small values it 
diverges and the processes which suppress its growth should be taken into account. 

But even for q/2p, > A when the integration should be performed over the whole velocity space 
it is not possible to obtain an analytical expression for the PTF in terms of elementary functions. 
However, we can calculate the coefficients W(q), W'(q), W"(q) and W'(q) expanding the integral 
into a series in the vicinity of the point (q-p) =0 & p 2 — when it is necessary. 

2 27-3 

np n b u 
W[q) 4p 2 qq 2 ^ 

X [7+ 1 q 2 1 + 2 
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where q = q/2/j,. 



W'(q) = 



Tip 



7 



27-3 
n b u ' 



4M 2 ii 2 q 7 + 3 ^2(7+2) 
A 2 3 2 



q 2 7 + 4 



7 + 2 



(49) 



W'(q) 



1 



2 27-3 
7T/T n b u ' 

Mm ^2(7+2) 

A 2 1 



W"{q) = 
'A 2 1 



7 + 2 
2tt P 2 



q 2 7 + 3 



7(7 + 1) n b u 



27-3 



M 2 ^ 2 g 7 + 3 q 2 ^+ 3 ) 



1 



(50) 



(51) 



g 2 7 + 4 7 + 2_ 

The convergences of all coefficients above at large values of q fulfills when 7 > 0. Remarkably for 
the dusty particles interacting with electrons the coefficient W'(q) is not equal to zero although it 
is small due to the mass ratio of electrons and ions. 

Similar analysis can be done for the dusty particles interacting with ions. 



4 Conclusion 

A more general Fokker-Planck-like equation is derived on the basis of the master equation. The 
equation is applicable for any values of the transferred momentum in a collision act unlike the usual 
Fokker-Planck equation which is valid for relatively small values of the transferred momentum. The 
mass ratio of the interacting particles is not important for the derived Fokker-Planck-like equation. 
It is based on the mere assumption that there are two species of the colliding particles, namely, 
with large and small velocity values. 

The coefficients for the general Fokker-Planck-like equation are calculated using the general 
expression for the probability transition function describing elastic collisions. The scattering particles 
were assumed to have a velocity distribution function with power-type tail. It leads to the power-type 
dependence of the PTF on the transferred momentum and the corresponding process of diffusion 
in velocity space has anomalous character. Three different examples of the differential cross-section 
have been considered. The hard-sphere model of collisions as well Coulomb model demonstrate 
similar behavior of the PTF. However, for Coulomb interactions PTF has a more strongly decreasing 
tail. The case of dusty particles scattering on plasma particles (electrons or ions) is more complicated 
leads to the less trivial dependence of the PTF on the transferred momentum. 
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